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Приводится методика расчета термограмм тепловой обработки вареных колбасных 
изделий с автоматической корректировкой параметров процесса. Даны начальные 
и конечные термограммы для колбас в натуральной и искусственной оболочках. 
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Введение. На сегодняшний день развитие производства практически лю- 
бой отрасли невозможно без применения информационных технологий 
(ИТ) — информационных взаимодействий, информационных структур и 
субъектов, осуществляющих сбор, формирование, распространение и ис- 
пользование информации. В не меньшей мере это касается и предприятий 
мясной промышленности. Современное предприятие, претендующее на 
выпуск конкурентоспособной продукции, должно не только оперативно ре- 
агировать на требования рынка, но также вести свою деятельность в соот- 
ветствии с все возрастающими требованиями по контролю качества и без- 
опасности мясных продуктов. 

В такой ситуации становится очевидно, что процесс контроля каче- 
ства не должен и не может быть ручным. Предприятия, на которых коэф- 
фициент автоматизации процесса контроля качества равен 100%, имеют 
колоссальное тактико-техническое, технологическое, качественное, а в ко- 
нечном итоге и экономическое превосходство над конкурентами, использу- 
ющими для тех же целей ручные ресурсы. 

В процессе производства мясных изделий наиболее критичным в 

спектре мероприятий, направленных на контроль качества, является этап 
тепловой обработки. Таким образом, ясно прослеживается необходимость 
мониторинга и прогнозирования характеристик продукции на указанной 
технологической операции. В настоящей статье приведена методика расче- 
та показателей процесса тепловой обработки на примере вареных колбас- 
ных изделий. Однако область применения можно расширить также и на 
мясные консервы. 
Методика эксперимента. При расчете термограмм тепловой обработки 
вареных колбасных изделий необходимо выделить две группы изделий: ва- 
реные колбасы в искусственных и натуральных оболочках. Принципиаль- 
ное различие между ними состоит в том, что их термограммы состоят из 
различного числа участков, соответствующих различным стадиям термиче- 
ской обработки. 

Для первых термограмма в общем случае состоит из шести 
участков: прогрев от начальной температуры до температуры варки; соб- 
ственно варка; снижение температуры до температуры охлаждения; вы- 
держка при этой температуре; снижение температуры до температуры до- 
охлаждения; выдержка при этой температуре [1]. 
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По существу, характер термограммы вареных колбас в искусствен- 
ной оболочке соответствует характеру термограммы мясных консервов, 
различие состоит только в продолжительности отдельных этапов термооб- 
работки и значениях температуры на каждом из них. 

Термограмма колбас в натуральных оболочках содержит десять 
участков: добавляются участки прогрева до температуры подсушки; соб- 
ственно подсушка; прогрев до температуры обжарки; собственно обжарка. 
Все последующие участки термограммы аналогичны термограмме колбас в 
искусственной оболочке. 

Автоматизированный метод расчета термограмм вареных колбас 
основан, на использовании математической модели процесса. В качестве 
физической модели колбасного изделия выбран цилиндр неограниченной 
длины. Ввиду симметрии формы цилиндра и симметрии нагрева рассмат- 
ривается сегмент в 2 цилиндра и полученный результат симметрично 
переносится на оставшиеся '2 цилиндра. 

Исходные данные: 

‚ расчетная температура и расчетное время прогрева централь- 

ной зоны изделия; 

‚ закон изменения температуры греющей среды; 

е Предельно допустимая температура поверхностного слоя изде- 
лия; 

‚ Геометрические размеры изделия; 

е Температура предварительного нагрева изделия; 

‚ физические и теплофизические параметры содержимого изде- 
лия (р-плотность; С - удельная теплоемкость; а - коэффициент 
температуропроводности; о- коэффициент теплоотдачи; „, - ко- 
эффициент теплопроводности). 

Принятые допущения: 

‚ задача решается в трехмерной постановке, но в силу осевой 
симметрии и значительного превышения длины цилиндра по от- 
ношению к его диаметру задача сводится к одномерной по ра- 
диальной пространственной координате; 

‚ принимается кусочно-линейный закон изменения температуры 
внешней греющей среды; 

. продукт представляется идеальной, непрерывной, однородной, 
изотропной средой с тонким, отличающимся по своим теплофи- 
зическим характеристикам внешним слоем; 

е  Температурное поле внутри термокамер симметрично относи- 
тельно цилиндра; 

‚ Геометрические размеры цилиндра в процессе термообработки 
не изменяются; 

‚ основные теплофизические характеристики продукта в общем 
случае зависят от температуры и изменяются в процессе термо- 
обработки. Однако экспериментальным путем выявлено, что ко- 
эффициенты а и Л для вареных колбас изменяются незначи- 
тельно (в пределах 4-5 %) и в расчетах могут приниматься за 
константы. 

Математическая модель. В силу того, что длина цилиндра значи- 

тельно превышает его диаметр, прогрев в основном происходит со стороны 
боковой поверхности, и тепловым потоком со стороны торцов можно пре- 
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небречь. Математическая модель, описывающая процесс распространения 
тепла в изделия цилиндрической формы, сводится к одномерному (по од- 
ной пространственной координате) уравнению теплопроводности Фурье в 
цилиндрических координатах. Из-за осевой симметрии нагрева можно 
рассматривать процесс распространения тепла только до центра батона. 

В качестве начального условия задается распределение температу- 
ры по толщине изделия в начальный момент времени. 

Технологическим операциям обработки колбасных изделий наибо- 
лее соответствуют граничные условия третьего рода, характеризующие за- 
кон конвективного теплообмена между поверхностью тела и окружающей 
средой. В этом случае количество тепла, передаваемого в единицу времени 
из окружающей среды через единицу боковой поверхности изделия внутрь 
его, прямо пропорционально разности температур между окружающей 
средой и поверхностью тела. 

Математической моделью, описывающей нестационарное поле тем- 
ператур, содержащих продукт, является параболическое уравнение тепло- 
проводности Фурье в цилиндрических координатах с переменными коэф- 
фициентами, с заданными начальными условиями (начальное поле темпе- 
ратур в продукте перед очередным этапом термообработки): 


2" А аа 





1 
91 С 92 г д (0 


0<г<В, (ЕВ ` 
Начальные условия (температура продукта до начала термообра- 
ботки): 
(("/! |= И, = сопб, 1=0, << В (2) 
Граничные условия (на боковой поверхности изделия): 
9 
| 5 Ив (Ё- Ц (2, | ‚ 9= — г=А, («КВ . (3) 
р 
Условия симметрии геометрии и нагрева в центре изделия: 
900: 
И 0. (о, | ‚ г=0, (ЕВ, (4) 
г 


где Ц - текущая температура продукта; Ц‹- начальная температура про- 
дукта; Ч» — температура греющей среды; Ц — конечная температура 
продукта; г - радиальная координата; В — внешний радиус изделия. 
Метод решения. Для решения задачи при изменяющейся во време- 
ни температуры греющей среды использовался метод конечных разностей 
(«неявная» схема). Возможность применения этой схемы обусловлена од- 
номерностью задачи по пространственной координате, при которой коэф- 
фициенты полученной системы линейных алгебраических уравнений обра- 
зуют ленточную, трехдиагональную матрицу, решаемую за небольшое чис- 
ло операций методом «прогонки». Выбор «неявной» схемы позволяет так- 
же существенно сократить время расчета за счет несвязанности шага по 
времени с шагом по пространственной координате. 
Выбрав шаг по пространственной координате Н и шаг по времени /, 
в результате замены дифференциальных операторов в уравнениях (1)-(4) 
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разностными операторами, приходим к системе линейных алгебраических 
уравнений: 


И0= Оо, 1=0 п - начальные условия; 











(5) 
а = 0 - условие симметрии; (6) 
Ч Ч а От 2 О я Он Чен 
| 12 (2-11 2 
1=1 п-1 г" 
Е. г Ве й В 
= р ЕН = 5 ИО ‚5 1-П, (8) 


где Ги п - индекс точки разбиения и число точек разбиения по про- 
странственной координате; К — индекс точки разбиения по оси време- 
Ни. 

Получив из уравнения (5) выражение для температуры в фиктив- 
ной законтурной точке (1=0) — Шок! и подставив его в первое из уравнений 
системы (6), приходим к следующей системе линейных (переход к линейно- 
сти возможен на основании допущения (6)) алгебраических уравнений: 

сады + Сы = Ь; 1=1, 


сна + Сына Сун Оныыа = В; 1=2 пв-1, (9) 
Ва С, „Она = Ь,; 1= п. 
Коэффициенты матриц С и В вычисляются по формулам: 
7 2а . 
т! р Яя 1= 1 
а[3- 21 12 а|[ 21-1 
а Е ее ЕЙ 1=2 п-1 
" 21 (1-1 И | 21 (1-1 (10) 
Е ЕП 
ь = -, ЕТ и БЕ ЖОЙ, т 


Полученная система линейных алгебраических уравнений (9) реша- 
ется методом «прогонки». 


(В 2) | 
От: ЕЕ т Он р 


) ) 


(7. =т+КО, 1=2:п+1 


ЬК+1 1 НР 
а = И, (11) 
где И 1=1 
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Е 1+1 
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сн. 1К;. |7 бы ь 


При прямом ходе метода «прогонки» рассчитываются коэффициен- 
ты т, К во всех точках разбиения по пространственной координате и 
определяется температура в краевой точке на границе продукт-греющая 
среда в К+1-й момент времени. При обратном ходе по рекуррентной фор- 
муле (см. формулу (11)) определяется температура во всех остальных точ- 
ках разбиения. По мере расчета поля температур по толщине батона кол- 
басы рассчитывается величина стерилизующего эффекта, который для ва- 
реных колбас проявляется при температурах выше 54°С. 
Начальные термограммы для вареных колбас в искусственных и на- 
туральных оболочках приведены в табл. 1, 2. 
Таблица 1 
Начальная термограмма для вареных колбас в искусственных оболочках 









































Температура, °С Тангенс 
= У(1)=15 {апг(1) = 0,08 
15г(3)=100 {апзг(3) = -0,05 
15г(4)=103 {апзг(4) = 0,00 
157 (5)=120 \(5)= 10 чапег(5) =-0,03 
15г(6)=121 апэг(6) = 0,00 
15г7)=150 {апзг(7) = 0,00 


Таблица 2 
Начальная термограмма для вареных колбас в натуральных оболочках 


вн: ВИ р ь 
МИН 


ЧапегС) = 0,750 
у(2)=60 {апзг(2) = 0,000 
{5"(3)=21 у(3)=60 {апзг(3) = 0,320 


Чапзг(а) = 0,000 
апзг(8) = 0,000 
{57(10)=127 {апзг(10) = 0,000 


{57(6)=180 м(6)= 5 чапзг(6) = 0,000 


Продолжительность этапов подсушки, обжарки, охлаждения и до- 
охлаждения определяется условием достижения требуемых температур в 
центре батона на каждом этапе термообработки (табл.3). Продолжитель- 
ность же этапа варки определяется из условия достижения требуемого сте- 
рилизующего эффекта. 





Таблица 3 
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Требуемые температуры в центре батона 


Наименование этапа _| Подсушка | Обжарка Охлаждение | Доохлаждение 





ра, °С 


Алгоритм расчета термограммы для вареных колбас в натуральных 
оболочках представлен на рис.1. Результаты решения представлены в 
табл.4—-7 и на рис.2,3. 
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Таблица 4 
Исходные данные по колбасам в натуральной оболочке 
































Коэффициент температуропроводности (а) 0,0015 
Коэффициент теплоотдачи (аНа) 0,004 
Коэффициент теплопроводности (1атЬада) 0,004 
Начальная температура (ч0) 15 
Продолжительность процесса (Тт) 14400 
Шаг разбиения по радиусу (п) 0,1 
Шаг разбиения по времени (1) 30 





Шаг печати по времени (респ) 
Радиус батона (г 
Число точек разбиения по ради 
Число участков термограммы (№ егт 
Нормативный стерилизующий эффект (Ю) 

























Таблица 5 
а расчетов для колбас в натуральной оболочке 


Е. ВЕН НИИ 
частка 


| Время, мин | 11| |278 28,5; | | 114 | 115 | 


Ба Я ый Вы Вы бы ыыы 
‚С 
В С И ОС С К С о ВС Е ОИ ВБ 


ЕЕНТИИЕ эффект Ёа=28,13501. 
























































Таблица 6 
Исходные данные по колбасам в искусственной оболочке 
Коэффициент температуропроводности (а) 0,0015 
Коэффициент теплоотдачи (аа) 0,004 
Коэффициент теплопроводности (1атБда) 0,004 
Начальная температура (ц0) 15 
Продолжительность процесса (Тт) 14400 
Шаг разбиения по радиусу (Н) 0,1 
Шаг разбиения по времени (1) 30 
Шаг печати по времени (респ) 10 
Радиус батона (г) 4,5 
Число точек разбиения по радиусу (п) 9 
Число участков термограммы (№егт) 7 
Нормативный стерилизующий эффект (о 18 
Требуемая температура в центре батона, °С 70 
Таблица 7 


Результаты расчетов для колбас в искусственной оболочке 


Вад ааа в а- 
ИИ ВИ В СЕВ ЗИ ИЕ ЕЕ ЕСИ 


я 15 
ра, °С 


[Тангенс Г озын [ о | 939 [0 | 98| | 
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Коэф. Температуропроводности а; Козф. Тепло- 
отдачи аа; Коэф. Теплопроводности |агада 
Радиус батона К; Число точек разбиения п 
Начальная температура 00; Конечное время ТМ, 
Шаг по времени 1: Шаг печати по времени оесь 






В=ааЛаша В=+-1) 
6[1,1]=-1; с[1,2]=1; с [уз-1]=5*5/2-1 
с[ву]=1+5* 5/2 5[1]=0 


с[1,1-1] :=а* (3-2*1) / (2**в* (1-1); 
с[1, 1] :=1/1+2*а/ в/в; 
СГ1, 1+11 :=а* (1-2%1)/ (2*Н*В* (1-1); 





К[ 1] :=-6[1,2]/[1,1]; 
и[1] :=6[1]/6[ 1,1]; 


КГ11 :=-сГ1,1+11 ИГ Г1,1-11*КГ1-11+сГ1,11) 


Расчет и печать 
площади под кривой 
конец 


9:=3+1; 
91:=91+1; 


Подпрограмма расчета 
температуры среды 


Рис.1. Алгоритм расчета термограммы вареных колбас 
в натуральных оболочках (окончание см. на с.117) 





117 


Вестник ДГТУ, 2009. 79. №1(40) 

















= А 


| ш(1]:=(6[1]-с[1,1-1]щ[1-1])/(с[1,1-1]*[1-1]4+с[41,1]); 


К[1] :=-с [1,141] /(6[1,1-1]*К [1-1] +°[1,1]) 


4[1]:=%1[1) 


Рис.1. Окончание 
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Рис.2. Термограмма вареных колбас в натуральной оболочке 
Стерилизующий эффект Ёа=18,27999 


Температура, С 


юр зря 


Время, мин 





— - -Тер- - - "Тов Тц 


Рис. 3. Термограмма вареных колбас в искусственной оболочке 


Полученные результаты послужили основой для разработки модуля 
автоматизированного расчета оптимальных параметров технологического 
процесса термообработки мясных изделий [2]. 
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